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Abstract 

We reconsider the possible presence of charge and colour breaking minima in the scalar 
potential of the minimal supersymmetric standard model (MSSM) and its minimal gener- 
alization with R-parity explicitly broken by bilinear terms (RMSSM). First we generalize 
some results previously derived for the MSSM case. Next we investigate how robust is 
the MSSM against its RMSSM extension. We examine the constraints on the RMSSM 
parameter space that follow from the required absence of charge breaking minima in the 
scalar potential. We point out the possibility of generating non-zero vacuum expecta- 
tion values for the charged Higgs field which is not present in the MSSM. However, given 
the smallness of neutrino masses indicated by neutrino oscillation data, we show that the 
RMSSM represents only a slight perturbation of the MSSM and is thus as safe (or unsafe) 
as the MSSM itself from unwanted minima in the scalar potential. 
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I. INTRODUCTION 



Softly broken supersymmetric models contain a fairly large number of scalar 
fields not present in the standard model. Their existence leads to a complicated 
scalar potential, which might contain undesirable minima which spontaneously break 
charge and/or color symmetry, a situation which can not happen within the Standard 
Model. The condition that the "realistic" minimum is the global minimum of the 
theory can be used to obtain restrictions on the parameter space of supersymmetric 
models, as already realized more than 20 years ago P, ^ • This way a disadvantage 
of supersymmetry may turn into a virtue by shedding some light into the unknown 
supersymmetry breaking mechanism itself. 

Due to the enormous complexity of the full scalar potential in the minimal su- 
persymmetric extension of the standard model (MSSM) early papers on this sub- 



ject PJ, 12, y, y] have only^nalyzed particular, but especially dangerous directions in 
field-space. Casas et al 5] have presented a more detailed analysis of this subject. 
They were able to show that in the constrained MSSM (CMSSM) with minimal 
supergravity boundary conditions strong constraints arise ruling out sizeable parts 
of the parameter space ^ . 

Similar studies in R-parity violating versions of the MSSM, however, have not 
been published ^. Our main goal is to present a detailed analysis of the 'unbounded- 
from-below' (UFB) as well as charge/colour breaking minima (CCB) in the bilinear 
R-parity breaking model (RMSSM) This model breaks lepton number and 
R-parity explicitly through the simplest bilinear terms. The justification for such 
emphasis is threefold. 

First, it represents the simplest possible scheme of R-parity violation, a mere 
six parameter extension of the MSSM. It is therefore interesting to investigate the 
"stability" of the MSSM against such "innocuous" perturbation. For this reason we 
can also call this model the generalized MSSM where R-parity breaks in the minimal 
way. Second, this model is motivated by the fact that it produces the paradigm 
for the idea that supersymmetry is the origin of neutrino mass leading to a 
pattern of neutrino masses |3] that successfully describes current neutrino data 1^ . 
Last, but not least, it represents the only model of R-parity breaking consistent 



^ The work of Abel and Savoy |g| contains a discussion on the possibihty of Ufting flat directions by 
adding cxphcit trihnear R-parity violating terms to the superpotential. However, they discuss 
the impact of bilinear terms only briefly. This is our main emphasis. 
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with a spontaneous violation of R-parity 0, 
fundamental theory, that breaks the symmetry. 



12 1, where it is the vacuum, not the 



In this model the atmospheric neutrino mass scale 131 is generated at t 



le tree- 



level, through the mixing of the three neutrinos with the neutralinos |1J], in an 
effective 'low-scale" variant of the seesaw mechanism. In contrast, the solar mass 
and ..ing= needed to accon,. fo. =o,a. neutnno data a.e generated .ad.a- 

tively Q. 

A very important difference between such a supersymmetric approach to the 
origin of neutrino mass and seesaw-type schemes, is that here the dimension-five 
operator responsible for (Majorana) neutrino masses is generated at an accessibly 
low energy scale - namely the weak scale. This makes this model potentially testable 
by experiment. In fact it has been shown that such a low-scale scheme for neutrino 
masses has the advantage of being testable also "outside" the realm of neutrino 
physics experiments. Although neutrino properties can not be predicted from first 
principles, interpreting current neutrino data in this framework imp lies unambigu- 
ous tests of the theory at accelerator experiments 0, Q, Q, boi 21 1 which can 
potentially be used to falsify the model. 

This paper is organized as follows. In the next section we will briefly recall 
some basics of the discussion on CCB and UFB bounds in the MSSM. This will 
serve as a basis for section 3, where we will discuss new features related to the 
R-parity violating terms. We show how the bounds from unbounded-from-below 
directions have to be modified, once non-zero bilinear R-parity violating (BRpV) 
terms are allowed. We point out the novel possibility to generate a non-zero vacuum 
expectation value of the charged Higgs field, albeit in regions of parameter space 
which are now excluded by neutrino physics We show that, given current 

data on neutrino masses, bilinear R-parity violation can be understood as a small 
perturbation of the MSSM. From the point of view of charge breaking minima the 
RMSSM is thus as safe (or unsafe) as the MSSM itself. We will then close with a 
short summary. 

II. REVIEW OF THE MSSM RESULTS ON UFB AND CCB 

To set up the notation, the superpotential of the MSSM can be written as 



W = Sab 



KQtU.H^ + hlQ\D,H^, + h%L\R,Hl - f^H^.H'^ . (1) 
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Here, /i^, /i^ and /i^ are 3x3 Yukawa matrices, Q, U and D are quark doublet 
and singlet superfields and L and R are the usual lepton doublet and singlet fields. 
Supersymmetry must be broken and the most general set of soft breaking terms 
allowed by the standard model gauge group under the assumption of lepton number 
conservation can be written as 

2 



+ 



i=l 

3 



-i ^ M,A.A. + Bat [j^i^QIU.Hl + A%h%Q\D,H^, + A%h%L\R,H^, 

i=l 

-B^iH^Hl) + /i.e. 
The Higgs doublets giving mass to the standard model fermions are 

and the parameters in Eq. ((21) are to be understood at some renormalization scale Q 
chosen to minimize the effects of the one loop corrections. This way we can neglect 
in the analysis the effect of the one loop radiative corrections j^. Without loss of 
generality, we now consider that the fields take the following vev's^, 

(if+>=0, {H^)=v., {H',)=v,, {Hl)=v^ (4) 

to obtain 

^Higgs = ("^L + ^^^) + ("^L + (^d + ~ ^BfiVuVd - ^g^vlvj 

+1 + g") K + + vt + 2vy_) + 1 - g'^) [vl + vl) vl (5) 

This Higgs potential has the property that f_ = 0. To see this we note that the 
potential can be written in the form, 

^Higgs = C^'"- + C^^- + C'o (6) 

where 

c. = \{g' + g") 

Co = I {g^ + g") {vl -vlY+ {m\^ + /i^) vl + {m\^ + /i^) vl - IB^v^v^ (7) 



(2) 



^ Our normalization here for the vev's differs from Refs. [lllSl by a factor of \f2. 



Now since g > g' we must have C2 > 0, unless + < ^. Therefore the 
minimum of the Higgs potential occurs for vanishing vev of the charged Higgs boson. 
By using the minimization equations, 

= -2Bfiv, + 2 {mj,^ + v^-]^ {g^ + g'^) {vj - vl) 

= -2B^v^ + 2 {ml^ + /i^) v, + ^ {g' + g'^) [v^ - vl) (8) 



one can find the value of the Higgs potential at the real minimum, 

Vmin = -U9' + 9"){vl-vlY (9) 



Eq. will be important to compare with the values of other (and potentially 
deeper) minima. 

Before starting the discussion of the dangerous directions, a word of caution 
should be added, namely, that the condition that the realistic minimum is the global 
one might actually be too conservative. In fact, it is possible that the universe 
resides in a false vacuum which is stable because the tunneling time into the global 
minimum is large with respect to the age of the universe. In this sense, CCB and 
UFB constraints on the supersymmetric parameter space are sufficient but might 



not be necessary, see for example 22, |23|. However, we will not follow this line of 
reasoning any further. 

A. UFB directions 

The 'unbounded-from-below' (UFB) directions are those where the quartic D- 
terms vanish and some coefficient (s) quadratic in the vev's are negative. Then the 
potential at the weak scale seems to be unbounded from below. However, this is a 
slight misnomer, since if one assumes that all soft masses are positive at the high 
unification scale, it appears that these dangerous directions are not really unbounded 
from below but there exists a true local minimum at some large scale. It then must 
be checked that this local minimum is not deeper than the physical one. As was 
shown in Ref. there are three kinds of such directions. The first and most obvious 
one corresponds to the D-flat direction where \vu\ = \vd\, all other vev's being zero. 



^ Casas et al. jsj assume that only + /i^ can be negative. Even though in mSugra at very 
large tan /3 values + /i^ < can occur in exceptional cases, we will follow their assumption. 
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The potential along this direction reads, 



VuFB-i = + + 2/x2 _ 2|i?^|) vl (10) 

and a sufficient condition to avoid developping a deep minimum at large values of 
the field is 

+ V-2|5/i| >0. (11) 

In principle, one should check the depth of the true minimum along the dangerous 
direction when this coefficient is negative. For simplicity, we will stick however to 
the condition given in Eq. (jllj) . 

The second dangerous direction corresponds to the case where a slepton Li takes a 
vev f j. Then a combination of Vu, Vd and Vi can cancel the D-term and the potential 
reads. 



which constrains the coefficient of the quadratic term as 



ml +fi^ + ml_ , *^^J ^>0. (13) 



Note that in the case of a universal mo at the unification scale the m/,. are usually 
the smallest soft masses at the weak scale. Dropping the universality assumption 
the bound obtained for m^,., Eq. (jl3|) . must be verified for the squark soft masses 
as well. 

Finally the last UFB direction corresponds to the case where Vd = but we have 
a neutral slepton Li with nonzero vev, like in the UFB-2 case. This direction is both 
D- and F-flat. The difference with respect to UFB-2 is that the F-term is canceled 
by giving vev's to the charged sleptons. The resulting potential reads 

VuFB-3 = {ml^ + ml) vl + f^ (ml + m\ + mM - -^-^ (14) 

Since m^^ must be negative in order to break electroweak symmetry and m|. is 
small when one assumes universality of the soft terms, the coefficient quadratic in 
Vu is generally negative. As shown in Refs. P, 0| in the case of universal soft masses 
at the GUT scale, the condition that the minimum along this UFB-3 direction is 
not deeper than the physical minimum implies mo > aMi/2, where a is a coefficient 
of C(l). 
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B. CCB minima 



For the classical CCB minima, dangerous negative contributions to the scalar poten- 
tial are generated by cubic {A-tjpe) soft supersymmetry breaking terms. Therefore 
these directions cannot be F-flat, but they are still D-flat. The traditional bound 
of Ref. jl| corresponds to the case where 

{Q') = {Hi) = {U) = v (15) 

all other vev's vanishing. This choice cancels the D-term and the potential reads, 

VccB = {3hy + 2AuKv + vn\^ + /x^ + m J + m^) (16) 

In order to avoid a very deep color and charge breaking minimum we must make sure 
that the parenthesis in Eq. ()16p never vanishes, which happens if the corresponding 
second order equation can not have real solutions. This leads to the well known 
condition, 

|A„p<3(m^^+/i2 + mJ + m2^) (17) 

A more complete and general analysis of this and similarly dangerous directions can 
be found in Ref. j^. Note again, that the bound given in Eq.lfTTj) for must be 
checked for all A-terms in the general non-universal MSSM. 

III. UFB AND CCB IN THE RMSSM 

The RMSSM is simply the bilinear R-parity violating model, defined by the 
following superpotential |7| 

W = WMssM + SabeiL^H'^ (18) 

and corresponding soft supersymmetry breaking terms, 

VsB = VMssA4 + B,e,L^Hl (19) 

It is therefore a rather mild extension of the MSSM. In the following it will be 
sufficient to consider for simplicity only a one generation version of the model ^. 



We do not believe that this simplification has any impact on the following discussion, since 
neutrino oscillation data require — <C 1 and intergenerational effects between different families 
of leptons due to BRpV terms scale as (^)^- 
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We are mainly interested in studying how the appearance of the new terms in the 
superpotential (and in Vsb) changes the conclusions which hold for the MSSM. Since 
the MSSM is the limit of the RMSSM when e — we expect that the results of the 
MSSM will hold in that limit. Note also that the structure of the trilinear terms is 
not modified, so conclusions like those of Eq. (fTTj) are expected also to hold in our 
case. Defining 



'X)=v,, {Hi) 



(20) 



one finds for the scalar potential 

V = Miy^ + Mjj^ [vl + vl) + Ml {v'^ + v'^) - 2BfiVdVu + 2B'eVuv' 



,9 (2 



12 



12 2 

V 



/2\2 



(21) 



where B' characterizes the soft supersymmetry and R-parity violating bilinear term. 
We note that it is not possible to have an UFB direction with non vanishing charged 
vev's in this potential, because the D-terms can not be made to vanish for f _ and 
v'_ different from zero. The minimization equations can be found in the usual way 
taking derivatives with respects to the fields 







12 



— {2eii — g^v^v'J) v — 1By,Vu 

9^ I 2 2 12 , 2 , l2\ , 9 / 2 2 I 

+2 (M|^ + /i^ + + 2 {B'ev' - B^iVd) 



12 



2 I /2\ y / 2 



'2 ^2 



Vd 



j2 



/2 2 

V — V_ 



9 I, .2 ..2 ,/2 , ,.2 ,./2\ 9 



n I T\/r2 , 2\ i/ / 2 2 /2 , 2 I2\ i/ / 2 

- (2e/i - g'^v-v'_) Vd + 2B'eVu 



.'2 



2/2 2 

Vd-v -v_ 



- [2efi - g^Vdv') v'_ 



9 /..2 ..2 I ,./2 I ..2 I ,./2\ 9 



o /' ;\ /r2 I 2\ I y / 2 2 , /2 , 2 , /2\ y / 2 2 i2 2 

- (2e/i - g\dv') v_ 



.'2\ 



(22) 



Since we are dealing with a set of five coupled equations this system is difficult to 
solve for the vev's. We can however use the following trick. Instead of solving for 
the five vev's we try to solve those equations for the three soft masses squared M'jj , 



12| and for the charged vev's. Using this approach we could find two 



Mfj^ and M| 
types of solutions. 

Before discussing the general case, however, we consider first the limit in which 
RMSSM is considered a perturbation of the MSSM. This is a reasonable approach 
since the BRpV parameters must be small to account for the neutrino data P|. 
Therefore we can pose the following question. Suppose that in the limit e — 
the parameters are such that the MSSM has no UFB directions or CCB minima. 
This means f n 7^ , f ^ 7^ and v' = V- = v'_ = If we now consider a small 
non-vanishing value for the e what will be the corresponding minimum? In order 
to answer this question in perturbation theory we write 

00 00 00 00 00 

vd = Y, ^\vu = Y, ^\ v' = Yl ^'^'^ ^- = ^- = 

i=0 1=0 1=0 2=0 i=0 

(23) 

Now we substitute back in the extremum Eq. ()22|1 and solve order by order in 
perturbation theory. The result that we get is as follows, 

(0) , (2) 2 , (4) 4 , 

Vd = vy +vye +vye +■■■ 

= + 42)^2 ^ ^(4)^4 ^ . . . 

y' = ^'(I)e + ^;'(3)e3 + ^;'(5)e5 + ... 



V- = Q 



(24) 



where Vu^ , v^^^ are the MSSM values for e = 0. This is precisely the solution of type 
I that we will discuss shortly. Note that if e 7^ then also v' 7^ 0. In fact, 

y' = tl^ g + (25) 

Ml-\{g^ + g'^) (t;f2_^(0)2j 

So we can formulate the following important result: // we start with the MSSM 
parameters such that in the limit e ^ the minimum has no UFB or CCB problems, 
then by turning on perturbatively a small value for e we get a correspondingly safe 
minimum of the RMSSM. However, as we will now discuss, in general there are two 
types of solutions for the minimum equations. 
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Type I 



This solution corresponds to the case where the charged vev's vanish. We are 
then in the situation studied usually in the bilinear R-parity model. We get 



efi fi +Bfi— + -V - v^) 



= -e -fi +Bfi Be [v^ - v - v^) 

]\j2 2 I f "^d o'^" A I fi'^ fi' ^ /' 2 /2 2\ 



v'_ = (26) 

This corresponds to the neutral Higgs potential that we will discuss further below. 
Here we just note that the value of the potential at the minimum can be shown to 
be 

VBn,v = -'-^{vl-vj-v-)\ (27) 



Type II 

In the general case we can find the solutions of the minimization equations in the 
following way. We start by solving the first three equations in Eq. (j^^ for the soft 
masses. We get, 

= MijO) - i + 9") {vl + V-) - 1 / ^ v'_v^ 

Ml = urn + 1 - - i + g'^) v'^_ -\g-^- v'.v_ (28) 

where M|f^(0), M^^(O) and M|(0) are the soft masses when f _ = f ^ = and are 
given in Eq. (j2El) • Now we substitute Eq. into the last two equations in Eq. 
to obtain, 

= -^2 I ^'2^_ _ v'vdv'_ + v'iv'_— - v.v'^ + 2eu (v^— - v'_ \ + 2Buv_— + g^v.vl 
V Vd J \ Vd J Vd 

r, 2 ( ,2 / / , 2 / ^' a\vd ^ ( v' ,\Vd , . 

y) = 9 f f _ — f VdV_ + v_v_ f _f _ 2eu f _ v_ — {29) 

\ Vd ) V \ Vd ) v' 
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Multiplying the second of the equations in Eq. (j29|l by v' jvd and adding them one 
obtains, 

v'_ = KV_ (30) 



where 

2B'e - g^v'Vu 

Finally we use Eq. (jHUj) to reduce either one of Eq. (pUj) to 



K = 31) 



= vl - Z^o) (32) 



where 



D2 = (k^ - —K 
Do = (v'^ - Vdv'K - vl[) - {Bvu + tv') — + 2e/iK (33) 

Vd 

Eq. (jH^ has the trivial solution t>_ = which corresponds to type I, the BRpV 
solutions. However, if 

^ > (34) 
we have a new type of solutions for the minimization equations. 




±\ —, v'_ = KV- (35) 
V D2 

As Do,2 do not have in general a well defined sign it can happen that such solutions 
do exist for some combination of the parameters. We will discuss this later in more 
detail. 

A. UFB Directions 

We have seen before that for the Higgs potential of the RMSSM the UFB direc- 
tions can only arise when the charged Higgs vev's vanish, otherwise it is not possible 
to cancel the quartic D-terms. The neutral Higgs potential obtained from Eq. (pT|l 
when t>_ = 0, f ^_ = is given by 



i2 

V 



2 /2 

-2Bii VdVu + 2B'e v^v' - 2fiev'vd + ^ {vl -vj- v'^) ' (36) 
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From this equation we can see that we can make the D-term vanish if we choose 
the condition 

vl = vl + (37) 
To implement this condition it is convenient to write 



Vd = Vu COS 6, v' = Vu sin 6 {31 



Then we get 



VNeutral = (39) 

where 

B{e) = \Mjj^ +e^ + fi^+ [Ml^ + /i^) cos^ ^ + (M| + e^) sin^ 6 

- 2B^ cos e + 2B'e sin 9 - 2^e sin 9 cos 9 (40) 

Therefore the condition for avoiding an UFB direction is that, 

B{9„,,,,) > (41) 

where ^min is the value of 9 that corresponds to the minimum of B{9). Now consider 
Eq. ()40p in the limit e ^ and take the derivative, 

^ = 2 sin ^ \-{Ml - MT) cos 9 + B^i] (42) 
d9 

The right hand side vanishes when ^ = and when cos^ = ,,-2 — ijT- These 
two solutions correspond to the UFB-1 and UFB-2 directions given in Eqs. flllj) and 
(Hn)), respectively, when e = 0. 

For e 7^ it does not seem possible to have an analytical expression for 9^i^. 
However for a given set of parameters it is always easy to verify whether Eq. (PT|) 
holds for 9 e [0, 27r]. It is also clear from Eq. (gOD that the MSSM condition, Eq. (ITT|l . 
is not enough to ensure that we are free from UFB directions. This fact can be best 
illustrated from figure (Q) that shows a typical example. 

One can see clearly that starting from a large value of -B(O) is not enough to 
decide upon the sign of i?(0min)- However it is easy to check numerically whether 
-B(^min) > or not. Therefore, although we lack a simple analytical formula, the 
criterium for avoiding UFB directions is easily implemented. 

Finally we comment briefly on the direction UFB-3. It can be easily shown that 
at large values of the field the potential in direction UFB-3 is given as 

VuFB-3 = {vfi\^ + + tB') vl^--- (43) 
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FIG. 1: B{0) as a function of 9 for an example where B{9mm) < but B{0) > 0. The 
right panel is an enlarged view of the left one close to the zeros of B{6). 



where the dots stand for irrelevant terms. Since in our notation eB' < this leads, 
in principle, to a slightly more stringent requirement than the one corresponding 
to the R-parity conserving MSSM. However, since ^ ~ 0(10"*^^""^^) is required by 
neutrino oscillation data j9|, this modification is_jiumerically irrelevant. This is in 
agreement with the argument presented in Ref. 



B. Nonzero charged Higgs and Slepton Vev's 

We now turn to the solutions of type II. We have already seen in Eqs. - 
that there are potentially dangerous solutions for the Higgs potential with nonzero 
vev's for the charged scalars. These solutions, if they exist, would provide new CCB 
solutions different from those already present in the MSSM, as explained above. 
As can be seen from Eq. fl34j) such solutions can exist if the parameters satisfy the 
relation Dq/D2 > 0, where the Di are given in Eq. ((221) • 

Since it does not seem possible to give a strict analytic criterion which relates 
the condition D0/D2 < (guaranteeing the absence of unwanted minima) to the 
parameters of the potential we have resorted to a numerical scan of the parameter 
space. Our approach to find the minima of the potential was as follows. We always 
started with a random set of parameters with zero charged vev's and subject to the 
requirement that, 

vl + Vd + = = (2v^Gf) = 174.1 GeV (44) 

Note that with this procedure we should always have, 

\v'\ 

1^1 = — < 1- (45) 

V 
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We then search for the global minimum numerically. If we find a minimum deeper 
than the realistic minimum but which breaks charge this part of parameter space 
should be discarded. Two examples are shown in Fig. (0). 



_ A' 
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10 
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FIG. 2: Range of RMSSM parameters where nonzero charged vev's for the Higgs and 
slepton fields are favoured over the reaUstic minimum for two examples of tan/3, left 
tan/? = 1.05, right tan/? = 1.2. Here we fix for convenience B = B' = fi = 100 GeV. For 
a discussion see text. 

The results shown in Fig. © can be understood qualitatively as follows. Starting 
with the definitions Eqs ()31|) and (jHS)) and taking into account the smallness of ^ 
one can show that in the limit e ^ we always have D2 > 0. On the other hand 
the condition Dq > requires 

(46) 



.0 9 tcLll 3 



1 ^ 25/itan2/? 



1 + tan^ /5 g tan /3 

Note that this condition is not strictly valid for tan/? = 1, because in this limit 
we can no longer neglect the terms proportional to e in the definitions of Dq and 
D2- Eq. (jmi) shows that charge breaking minima in the limit of small values of e 
require that v' take up a sizeable fraction of v. This trend is clearly visible from 
Fig. 0. The figure also illustrates how these solutions disappear very quickly with 
tan/? greater than 1. 

Although we find it amusing that such solutions exist, we wish to stress that 
consistency with neutrino data requires ^ ~ 0(10^*^^"^^) and ^ ~ (9(10^*^^^^^). 
We therefore conclude that the RMSSM is automatically safe from these unwanted 
minima in those "physical" parts of parameter space which account for the neutrino 
oscillation data. 
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IV. CONCLUSIONS 



We have studied charge breaking minima and unbounded from below directions 
within bihnear R-parity breaking supersymmetry. Such a "reference model" is 
nothing but the simplest broken R-parity version of the Minimal Supersymmet- 
ric Standard Model. We have first generalized some results obtained previously in 
the R-parity conserving MSSM. Subsequently we discussed new ways to generate 
a nonzero vacuum expectation value of the charged Higgs and slepton fields. How- 
ever, such unwanted solutions occur only in regions of parameter space which are 
now excluded by neutrino oscillation data. 

In summary it can be said that, given the data on neutrino masses, bilinear R- 
parity violation can be understood as a small perturbation of the MSSM. From the 
point of view of CCB and UFB directions the RMSSM is as robust as the R-parity- 
conserving MSSM: it is equally safe from unwanted minima in the same portions of 
parameter space. 
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